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The lattice summations of the potential energy of importance in the graphite system have been computed 
by direct summation assuming a Lennard-Jones 6-12 potential between carbon atoms. From these summa- 
tions, potential energy curves were constructed for interactions between a carbon atom and a graphite 
monolayer, between a carbon atom and a graphite surface, between a graphite monolayer and a semi- 
infinite graphite crystal and between two graphite semi-infinite crystals. Using these curves, the equilibrium 
distance between two isolated physically interacting carbon atoms was found to be 2.70 a, where a is the 
carbon-carbon distance in a graphite sheet. The distance between a surface plane and the rest of the crystal 
was found to be 1.7% greater than the interlayer spacing. Theoretical values of the energy of cohesion and 
the compressibility were calculated from the potential curve for the interaction between two semi-infinite ^ 
crystals. They were A£c=—330 ergs/cm* * and |8 =3.18X10“^ cmVdyne, respectively. These compared 
favorably with the experimental values of A£c= —260 ergs/cm* and /3= 2.97X10"^ cmVdyne. 


P HYSICAL properties of solids such as energy of 
cohesion and compressibility can be computed from 
appropriate summations of the potential energy inter- 
actions of the atoms in the crystal. Graphite is a par- 
ticularly suitable subject for such computations be- 
cause of its crystal structure. The large spacing between 
graphite layers greatly simplifies the calculations and 
makes it necessary to consider only the interlayer 
cohesion. Also, the same set of lattice summations can 
be used in calculating the compressibility and the 
energy of cohesion. Furthermore, exp>erimental values 
of both compressibility and energy of cohesion are 
available so that the results of the calculations can be 
compared to experiment. 

Several studies have been made involving the sum- 
mation of the potential energy function of carbon over 
the graphite lattice. Barrer,^ and Crowell and Young^ 
have used these summations to compute heats of ad- 
sorption on a graphite surface. Kraus® computed the 
energy of adhesion of carbon tetrachloride to graphite 
assuming that each graphite layer could be treated as 
having a uniform density and then summing over the 
layers. Brennan^ obtained a value for the interlayer 
binding energy of graphite from a quantum mechanical 
calculation of the repulsive energy between layers. He 
also assumed that each layer could be treated as a 
continuum but that summation over the layers was 
necessary. 

It is the purpose of this work to evaluate the lattice 
sums on which all studies of the type mentioned above 
must be based, to determine the potential energy func- 
tions important in the graphite system, and to apply 
these results to a calculation of the energy of cohesion 
and compressibility of graphite. 

^ R. M. Barrer, Proc. Roy. Soc. (London) A161, 476 (1937). 

* A. D. Crowell and D. M. Young, Trans. Faraday Soc. 49, 1080 
(1953). 

*G. Kraus, Progress Report No. 2, University of Cincinnati, 
Army Air Force Contract AF33(616) 231, 1953. 

* R. 0. Brennan, J. Chem. Phys. 20, 40 (1952). 


POTENTIAL ENERGY FUNCTIONS IN TERMS 
OF LATTICE SUMMATIONS 

It will be assumed that the potential energy between 
two physically interacting carbon atoms in graphite is 
given by the Lennard-Jones expression 


where r is the distance between carbon atoms and A 
and fo are constants. For the purposes of this paper it 
is convenient to measure distance in terms of a, the 
bond distance between carbon atoms in a graphite sheet. 
Then, letting y=r/a, Eq. (1) becomes 



In the graphite structure, two types of equilibrium 
position are possible for a carbon atom above a graphite 
sheet: one in which the^tom is above the center of a 
hexagon, and the other in which the atom is above the 
corner of a hexagon, i.e., another atom. In forming 
lattice sums, therfore, two types of summation are 
necessary. 

In the following discussion, an atom located above 
the center of a hexagon of a graphite sheet will be called 
an h atom and an atom located above the corner of a 
hexagon will be called a c atom. 

Let = be the summation of 1/y”* 

i 

terms where x is the distance between an atom above 
the center of a hexagon and the graphite sheet and y,- 
is the distance between the atom above the graphite 
sheet and the ith. atom in the graphite sheet. 

The potential energy function between an h atom 
and a graphite monolayer sheet is given in terms of the 
above summation as 

*(*) = -li (yo)«5, -S, *«(x)]. (3) 

a* 
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Table I. The lattice sums for the graphite system. 


X 


Sl,c^(x) 

Seo.k^(x) 

5oo.c®(a7) 

5(0,03® (x) 


Sl.c«(x) 


0.5 

3.19216 

65.78386 

3.212517 

65.80415 

69.06384 

1.57397 

4096.79209 

4098.366 

1.0 

0.83067 

1.519795 

0.8417850 

1.53019 

2.399910 

0.094182 

1.048563 

1.142745 

1.5 

0.220744 

0.2483125 

0.22740 

0.254972 

0.499582 

0.0052127 

0.0106040 

0.0158167 

2.0 

0.074605 

0.076071 

0.078876 

0.080342 

0.170715 

0.00041437 

0.00050398 

0.00091835 

2.5 

0.030890 

0.030987 

0.033780 

0.033877 

0.075836 

0.490619X10“-* 

0.515632X1(F^ 

0.000100625 

3.0 

0.014923 

0.014928742 

0.016963 

0.0169687 

0.039992 

0.81281 XIO-*' 

0.822298X10“** 

0.1635108X10-^ 

3.5 

0.0080577 

0.00805629 

0.0095477 

0.0100917 

0.0242800 

0.17505 XlO-s 

0.175493X10“** 

0.350543 X10“5 

4.0 

0.0047235 

0.00472194 

0.0058429 

0.0058429 

0.015158 

0.4611 X10“® 

0.46133 XKT** 

0.92243 XlO-6 

5.0 

0.0019348 

0.00193324 

0.0026111 

0.0026111 

0.0074934 

0.495 XlO-7 

0.44496 XlO-7 

0.93996 XIO-T 


The potential energy between an h atom and a semi- 
infinite graphite lattice is obtained by summing Eq. (3) 
over all the .graphite layers in the semi-infinite crystal. 
The result is 

(4) 

where 

i=o . • 

and 

*S*oo, = (6) 

t=0 

I is the distance, between monolayer planes in the 
graphite lattice. 

An analogous set of equations hold for interactions 
involving a c atom identical in form to Eqs. (3), (4), (5), 
and (6) except that the . subscripts c and h . are inter- 
changed. 

The potential energy function for interaction between 
a graphite monolayer sheet and a semi-infinite graphite 
crystal ^oo.«(^) is therefore given by 

i ■ . 

C^oo.c(^)+^oo,fcW],- (7) 

2(T 

where cr is the area, occupied by a carbon atom in a 
monolayer sheet so that 1/2 (t is the number of atoms 
of each type per square centimeter in the monolayer 
and ^ 00 , « is given per cm^ of surface. 

Finally, the potential energy function for interaction 
between two semi-infinite graphite crystals is obtained 
by adding Eq. (4) and its analog for a c atom (i.e., 
with interchanged subscripts) multiplying . by the 


Table II. CrowelFs values of the lattice summations. 


Values obtained 
by direct 
summation - 
.X ■■ 5oo.a® 

5«,c« 

5oo,a®« by 
approximation 
method 

(This 

work) 

5«.c® 

(This 

work) 

1.5. 


' 0.2540 

0.246 

0.2274 

0.2550 

2.0 

0.07773 

0.07971 

0.0799 

0.07888 

0.08034 

2.5 

0.03330 

0.03343 

0.0336 

0.03378 

0.03388 

3.0 

0.01620 

0.01666 

0.0170 

0.01696 

0.01697 

3i5i 


0.009283 

0.00959 

0.00955 

0.01009 


number of atoms of each type {h and c) per cm^ of 
graphite monolayer, and summing over all the layers of 
the second semi-infinite crystal. The result is 

A 

Lh(yoys^.«^^{x)-s^j{x)2 ( 8 ) 

2(to« 

where 

S^„^{x) = f: “(*+//)] (9) 

i=0 

and 

5^»®W = i:C.S„.*'>(s;+t7)+5.,/(:r+iO]. (10) 

t=0 

EVALUATION OF THE LATTICE SUMMATIONS 

The sums Si,h^ and were computed by direct 
summation for all atoms within a distance of 30a from 
the atom above the graphite plane for values of x 
ranging from 0.5 to 5. The contribution of the remain- 
ing atoms in the plane was obtained by integration and 
affected only the sixth decimal place of the summations. 
The sums Si,h^^ and Si,c^ were computed by direct 
summation for all atoms within a distance of 10a from 
the atom above the plane. The contribution of the 
remaining atoms in this case amounted to only 10“^^ 
and was therefore neglected. 

The sums 5oo,c®, and 5oo,oo® were obtained by 

carrying the summation over the sheets from i = 0 to 30. 

Because of the rapid fall of l/y^^ with distance it was 
found necessary to use only summations over adjacent 
layers to obtain *Sco,c^^ and 5oo,oo^^ 

In the application of Eqs. (6), and (7) to the calcula- 
tion of the other sums, Si, and ^i, c® were needed at 
values of x for which direct summations had not been 
carried out. For values of x below 3.5, values of the 
lattice sums were calculated from Newton’s interpola- 
tion formula. The interpolation was applied to log5i,/i® 
and log5i,c® rather than to Si,h^ and 5i.c® because the 
log plots of these functions are nearly linear and can 
therefore be interpolated more accurately than the 
functions themselves. For values of x greater than 
3.5, Si, and Si, c® were computed by integration over 
the graphite plane assuming a continuous, uniform 
density. 

For x greater than 12 it was sufficient to compute 
5oo, A® and Soo, c® for use in Eq. (9) by integration assum- 
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ing a continuous, uniform, three-dimensional density. 
The errors introduced by these procedures were 
negligible. 

The results of these calculations are listed in Table I. 

An approximate method of obtaining lattice sums for 
graphite has been given by Crowell.^ He evaluated the 
lattice sums of 1/y® for the interaction between an 
atom and a semi-infinite graphite lattice by the ap- 
proximate method and compared them to the results 
of direct summation for several positions of the atom 
above the graphite surface. His results lead to the 
values of a® and S^o, c® given in Table II. A compari- 
son of Tables I and II shows that Crowell’s values differ 
somewhat from those presented in this work. In every 
case Crowell’s values are lower than the corresponding 
summations obtained in this study. The source of this 
discrepancy is not immediately evident and can be 
revealed only by a step-by-step comparison of the two 



Fig. 1. The potential energy function for interaction 
between two carbon atoms. 


sets of calculation. The fact that Crowell summed over 
only the nearest 100 carbon atoms can account in part 
for his lower results but probably not completely. 

The lattice sums obtained by the approximate 
method agree equally well with both sets of direct 
summations. 

POTENTIAL ENERGY FUNCTIONS 

Before the lattice sums can be applied to the calcula- 
tion of measurable quantities, the values of A and yo 
in Eq. (2) must be known. Kraus^ has computed a 
value of ^ = 24.3X10~®® erg cm® by Kirkwood’s 
formula using constants given by Barrer,^ and yo may 
be obtained from the condition that at equilibrium, 



® A. D. Crowell, J. Chem. Phys. 22, 1397 (1954). 



Fig. 2. The potential energy function for interaction between a 
graphite sheet and a semi-infinite crystal. 


At equilibrium x=l=2M and differentiation of Eq. 
(10) gives 


/^‘^oo,oo®\ / /dSoo, oo^^\ 

yo®=2( ) / ( ) . (11) 

\ dx / 2.36' V dx / 2.36 


Calculation from Eq. (11) gave yo= 2.70^, The deriva'r 
tives were obtained graphically from plots of log5oo,«,® 
and log5oo,oo^^ y. Several values of the lattice 
sums between x=2 and x=3 were calculated by New- 
ton’s interpolation formula. <p, h, <pl. c, <Po,. a, : <Poo, ^ 
^ 00 , s, and ^ 00,00 were computed using the following set 
of constants : 


^ = 24.3X10“®® erg cm® yo=2.70 
a= 1.42X 10~® cm <r= 2.62X 10~^®,cm^ ; ... 

The potential energy curves for and ^oo.« ^re 

shown in Figs. 1 to 3. The equilibrium distance for each 
potential function was obtained from an expanded plot 
in the region of the minimum. These distances, given 
in Table III, require some comment. The equilibrium 
distance of approach between two carbon atoms is 2.70. 
This, of course, is for a hypothetical case; namely, two 
isolated carbon atoms that are interacting physically 
and not chemically, and for which the potential energy 
constants are the same as those for graphite. If a carbon 


Fig. 3. The poten- 
tial energy function 
for interaction be- 
tween two semi- 
infinite graphite 
crystals. 
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Table TII. The position of the minimum for the various poten- 
tial energy functions — distance measured in units of the carbon — 
carbon distance in a graphite monolayer. 


Potential function 

Equilibrium distance 

<p{x) 

2.70 


2.41 

<Pl.h(x) 

2.37 

^«o,c(^) 

2.40 

*P^h{x) 

2.37 


2.40 

^«0k»(^) 

2.36 


atom interacts with a graphite monolayer, the equi- 
librium distance decreases to 2.41 for a c position and to 
2.37 for an h position. This is as expected because of the 
increased number of atoms interacting with the atom 
above the sheet. Similarly, one might expect a further 
contraction in equilibrium distance for a surface layer 
on a semi-infinite crystal because of the interactions 
between the surface atoms and those in the layers 
beneath them. This interaction has/ been found by 
Lennard-Jones and Dent® to result in a 5% contraction 
in the surface layer spacing of alkali halides. However, 
the large interlayer spacing characteristic of the graphite 
structure leads to an expansion of 1.7%. The large inter- 
layer spacing results in a negligible contribution to the 
repulsive potential for all except the adjacent layers. 
Therefore, the repulsive potential between a plane and 
a semi-infinite crystal is essentially the -same as the 
repulsive potential between two semi-infinite crystals. 
However, the attractive potential does not decrease as 
rapidly with distance so that planes at larger distances 
also make an appreciable contribution. Thus, the at- 
tractive potential is less for interaction between a 
surface plane and a semi-infinite crystal than for inter- 
action between two semi-infinite crystals; The net 
result is an increase in interplanar spacing for surface 
layers. 

COMPARISON WITH EXPERIMENTAL DATA: ENERGY 
OF COHESION AND COMPRESSIBILITY 


1027 ergs/cm^ so that our results are in serious dis- 
agreement with Brennan’s regardless of which wave 
function is chosen as correct. Recently, experimental 
energies of cohesion have become available^ from heat 
of wetting ‘ data. The mean of eight such values for 
graphite gives A£c=— 260 ergs/cm^ with an average 
deviation from the mean of about 11%. A comparison 
of these values is summarized in Table IV. The agree- 
ment between our value and the experimental value is 
quite satisfactoiy considering the experimental error 
and the possible error in the value of A. Furthermore, 
the lattice sums yield a value of A£c that is valid only 
at 0°K whereas the experimental value was determined 
at room temperature. Since the magnitude of A£c de- 
creases as the ternperature increases, the discrepancy 
between the experimental and theoretical values is 
easily accounted for. 

The compressibility /3 can be computed from the 
relation 


1 d^U 
-=F— , 
dV^ 


( 12 ) 


where V is volume and U = energy of the crystal. If it is 
assumed that the total volume change on compression 


Table IV. Comparison of various values of the energy of cohesion. 


Method 

Reference 

AEc ergs/cm* 

Lattice summation 

This work 

-330.0 

Zener’s wave function 

4 

-69.5 

Kohlrausch’s wave function 

4 

-1027.0 

Experimental (heat of wetting) 

7 

-260.0 


is the result of changing the interlayer spacing and that 
the entire energy change is the change in the lattice 
energy ^oo.oo, then Eq. (12) becomes 

dl^ ' 


The energy of cohesion for graphite is given by the 
interlayer binding energy and may be calculated from 
the potential energy curves. If there were no difference 
in lattice spacing between a surface layer and interior 
layers, the energy of cohesion would be given by v?oo,oo 
(2.36) i.e., the value of the interlayer potential at 
equilibrium. This value, however, requires a small cor- 
rection because of the expansion of the surface layer 
spacing from 2.36 to 2.40. The energy of cohesion AEc 
is then given by 

AE,= ^oo.oo(2.36)-[^«..(2.40)-^^..(2.36)] 

therefore AEc= —330 ergs/cm^. 

Brennan^ has given two values of AEc based on two 
different wave functions. These are —69.5 ergs/cm^ and 


where /= interlayer spacing. The derivative {dJ^(p^,^/dP) 
was computed graphically from plots of log5oo,oo® and 
log5«,.oo^^ as a function of x. The result thus obtained 
was 3.18X 10“^^ cmV dyne. This is in excellent agreement 
with the experimental value 2.97 X10“^^ cmVdyne that 
Brennan® calculated by extrapolating Bridgeman’s® 
data to zero pressure. 

SUMMARY 

The various lattice sums of importance in the graphite 
system were computed and used to calculate the poten- 
tial energy functions for interactions between two 
carbon atoms, between a carbon atom and a monolayer 
graphite sheet, between a carbon atom and a semi- 
infinite crystal, and between two semi-infinite crystals. 


® J. E. Lennard-Jones and B. M. Dent, Proc. Roy. Soc. (London) 
121, 247 (1928). 


’L. A. Girifalco, Ph.D. thesis, University of Cincinnati (1954). 
8 P. W. Bridgeman, Proc. Am. Acad. Arts Sci. 76, 9 (1945). 
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